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The purpose of this work is twofold: to present mathematical expressions for the kinematics of 
an articulated mechanism and to perform numerical experiments with the implemented code. The 
system of rigid parts is made of two slender bars and a disk. In the original configuration, a constant 
counterclockwise rotation rate is imposed on the disk. In the modified version, this angnlar velocity 
varies linearly with the rotation angle to produce an average rate that is nearly the same as the 
constant case. Angles, velocities and accelerations are analyzed for a 90° tnrn of the disk. The 
numerical solutions show the inversion of the linking bar sense of rotation along with the start of 
deceleration for both bars. The paper and pencil solution of the original problem that may lead 
to a wrong conclusion is explained. Equations are derived from first principles and the code is 
placed under version control. Those in charge of vector dynamics courses may Hnd it useful as a 
project-based learning activity. 


I. INTRODUCTION 

The rigid body is an important idealization of moving and interacting parts and beings. The key feature is the 
fixed distance between two points in the body, since forces and torques cannot cause deformation of its shape [1]. 
The subject is relevant for the development of video-game physics engines since the simulation of motions and 
interactions of rigid bodies approximates reality fairly well [5]. Prosthetic limbs and robotic devices are other 
important application areas of this theme. 

Here we offer a tool to stimulate engineering students to continue their learning path in the subject. We show 
mathematical formulations and numerical simulations of the interdependent motions. 


II. MATHEMATICAL FORMULATIONS 

The analysis has its origins in [3], where there is a solved example for the instantaneous position as shown 
in Fig. ^ and general guidelines for a computer project are proposed. In what follows, equations for the kine¬ 
matics are derived showing angle and length constraints, as well as velocity and acceleration relations while the 
disk undergoes either a constant 2 rad/s, or a variable (0.10 -I- 2.3 • a) radjs, with a being the turn angle of the 
disk [Fig. [1(b)] . 


A. Angles and length constraints 

The three angles describing rotation of the articulated mechanism are a, j3, and 7. A one letter notation is 
introduced to write trigonometric functions as: 


1 = 

sin a 

m 

= cos a 

(1) 

n = 

sin p 

0 

= cos p 

(2) 

P = 

5m 7 

Q 

= cos 7 

(3) 


Since the length of the linking bar AB cannot change, the horizontal displacement of pin A must be equal to the 
horizontal displacement of pin B, or: 







(a)Initial position 



(b)Intermediate position 


FIG. 1: Articulated mechanism at its (a) initial and (b) intermediate positions. Reproduced from [3]. 


OA - I = BC - BC • o 


(4) 


An overbar means the length of a given straight line segment. It can be inferred from Fig. 1(a) that OA = BC. 
Simplifying: 


1 = 1 - 


(5) 


A right triangle with AB as the hypothenuse will have, at its intermediate position in Fig. 1(b) the following 
opposite side to angle 7 : 


BC■ n + OA■ m 


( 6 ) 


And the adjacent side: 


OC - OA - I - BC- 


(7) 


Applying the Pitagorean theorem in Fig. 1(a) 


AB = ^(OC - BC)^ + {OAy 


( 8 ) 


These relations are needed to obtain /3, 7 , and the equations for velocities and accelerations in the following 
sections. 


B. Velocities 

With pin A viewed as belonging to the disk, its velocity is: 

Pa = w X poA (9) 

Where uj = 2 k rad/s for the constant rotation case, or w = ( 0.10 + 2.3 • a) k rad/s for the variable rotation 
one. For a going from 0° to 90° counterclockwise starting at the lowest position, w varies from 0.10 rad/s to 3.7 rad/s 
when Lo varies linearly with a. This accounts for a 1.9 rad/s average angular velocity. We can thus write: 




















a{t) = 1.9 • t 


( 10 ) 


Taking the time derivative of uj we get: 


uj = 2.3 • a 


( 11 ) 


And 


w = (2.3) (1.9) = 4.4rad/s^ 
With pin B belonging to the output bar BC, its velocity is: 


( 12 ) 


Vb = ^BC X PCB 


(13) 


Where ixibc is the still unknown angular velocity vector of bar BC. The motion of linking bar AB can be described 
using Chasles theorem [5]: 


Va = Vb + l^AB X Pba (14) 

Where t^AB is the still to be determined angular velocity vector. Va and Vb are always tangential to the circular 
trajectories of pins A and B. For the mechanism’s intermediate position shown in Fig. |l(by| x and y components for 
both velocities can be written as: 


VAcc = m- Fa VAy = l- Fa (15) 

Vbx = n - Vb VBy = o - Vb (16) 

The fixed vectors in the rigid bodies going from O to A (poa), C to B (pcb), and B to A (pba) can also be 
expressed using their x and y components. An inspection of Fig. |l(b)| allows one to write: 


POAx — ■ OA 

POAy = -m ■ OA 

(17) 

PCBx = —0 ■ BC 

PCBy = • BC 

(18) 

PbAx = -q- AB 

PBAy = -p- AB 

(19) 


Substituting 


(15) - (19) back in Eqs. ® , (131, and (141: 


and 


Vb = Fa 


m ■ q + I ■ p 
n • q + o ■ p 


( 20 ) 


^AB — 


AB 


m ■ o — I ■ n 
o • p + n ■ q 


( 21 ) 


With (201 in Eqs. (131 and (16), wbc can be calculated. 















C. Accelerations 


We now need to distinguish between the constant and the variable uj scenarios. Pin A only has radial acceleration 
with constant w, but it has both radial and tangential components otherwise. Pin B, as before, will have both radial 
and tangential components. 


1. Constanta) 


Taking the derivative of Eq. (14) with respect to time: 


a A = Sbr + asT + <^ab x pba + ^^ab x {wab x pba) 

The linear accelerations aA^ Sbr and qbt can be expressed in terms of unit vectors i and j: 


( 22 ) 


aA = w 


OAj ■ (^—l i + mj^ 


(23) 


Sbr = {<^BC^ BC) ■ (^oi - nj'j 


(24) 


Sbt = - (ujbcBC) ■ (^ni + oj'j 


(25) 


Substituting back into Eq. (22) and working the vector algebra, one is able to obtain expressions for angular 


accelerations of bars AB and BC. As these expressions are long, they are grouped into four different terms as ABi 
through AB 4 for cvab, and BCi through BC 4 for cvbc- It should be noted that these quantities point along the 
z-direction. The final result is: 


ABi 

= OA ■ {lo + mn) 

AB2 

= u)BC^ BC-{o^+n^) 

( 26 ) 

AB3 

= loab"^ AB ■ {oq - np) 

AB4 

= AB ■ {op + nq) 

( 27 ) 


bjAB = 

ABi + AB2 + AB:^ 


( 28 ) 


AB4 


BCi 

= uP' OA ■ (mp — Iq) 

BC2 

= ujab'^ AB ■ (jp^ + q^) 

( 29 ) 

BC3 

= lobc^ BC ■ {np - oq) 

BC4 

= BC ■ {nq + op) 

(30) 


BCi - BC 2 + BC 3 
1^4 


UJbC = — 


(31) 






2. Variable uj 


Since the disk now has a variable rate of rotation w, pin A’s acceleration is: 


a A = Sar + a AT 

Referring back to Fig. |l(b)[ the tangential component can be written as: 


(32) 


Sat = w OA 


{mi + 


(33) 


The radial component is still given by Eq. (23). As for the constant case, angular accelerations for bars AB and 
BC are written as: 


ABi = up' OA ■ {lo + mn) 
AB 3 = ujbc^ BC ■ (o^ + n^) 
AB^ = AB ■ [op + nq) 


AB 2 = u! OA ■ {mo — In) (34) 

AB 4 = ujab'^ AB ■ {oq — np) (35) 

(36) 

(37) 


Cab 


—ABi + AB2 — AB^ — AB4 

' AB 5 ■ 


And, 


BCi = uP OA ■ {mp — Iq) 
BC 3 = UJBC^ BC ■ {np - oq) 
BC 5 = BC ■ {nq + op) 


BC 2 = up OA ■ {mq + Ip) 
BC 4 = UJAB^ AB ■ 


(38) 


(39) 

(40) 

(41) 


Cbc 


-BCi - BC 2 - BC 3 + BC 4 


(42) 


This concludes the kinematical formulation of the articulated mechanism that has been implemented numerically 
and placed under git version control [ 1 ]. 


III. NUMERICAL EXPERIMENTS 

A. Constant uj 

Geometrical and kinematics parameters for the 90° turn of the disk are shown in Fig. Of note is the change 
of rotation sense for the linking bar AB from counterclockwise (positive ujab) to clockwise (negative ujab) when 
a = 30° (top right panel). This feature is also evident when 7 reaches a maximum (top left panel), and Vb and 
uiBC remain nearly constant around a = 30° (bottom left and top right panels). Angular accelerations, despite their 
intricate mathematical expressions [Eqs. ( [2^ through are slowly varying functions after the start-up period. 









a (degrees) 



FIG. 2: Top left', p and 7 vs. a. Bottom left-, linear velocities of pins A and B vs. a. Top right: angular velocities of bars AB 
and BC vs. a. Bottom right: angular accelerations of bars AB and BC vs. a. Angular step size Aq = 0.5°. 

TABLE I: Head and tail of linear and angular velocities. See text for comments on cells with asterisks. 


an 

pn 

7r) 

Va (m/s) 

Vb (m/s) 

u) (rad/s) 

UAB (rad/s) 

OJBC (rad/s) 

0 

0 

16.7 

0.6 

2* 

2* 

2* 

-6.7 









90 

90 

16.7 

0.6* 

0.18 

2 

-0.6* 

-0.6* 


Table |TT| has the computed angular accelerations’ initial and final values. The results for a = 0° serve to validate 
the numerical implementation since the computations match the textbook calculations in [5]. The same is true for 
the angular velocities of Tab. |T] 


TABLE II: Head and tail of angular accelerations. 


«n 

P(1 

7 0 

dJAB (rad/s^) 

ojbc (rad/s^) 

0 

0 

16.7 

-57.8 

192.6 






90 

90 

16.7 

0.0 

5.2 

























































B. Variable uj 


Results are now displayed in Fig. Note the linear increase of Va (bottom left panel), and how angular speeds 
and accelerations evolve without resemblance to the constant w case. In particular, it is now easy to see when both 
bars start decelerating during the 90° counterclockwise turn. 






FIG. 3: Top left: /3 and 7 vs. a. Bottom left: linear velocities of pins A and B vs. a. Top right: angular velocities of bars AB 
and BC vs. a. Bottom right: angular accelerations of bars AB and BC vs. a. Angular step size Aa = 0.5°. 


TABLE III: Head and tail of linear and angular velocities. See text for comments on cells with asterisks. 


o(°) 

pn 

7 n 

Va (m/s) 

Vb (m/s) 

u) (rad/s) 

LJAB (rad/s) 

wbc (rad/s) 

0 

0 

16.7 

0.03 

0 . 1 * 

0 . 1 * 

0 . 1 * 

-0.33 









90 

90 

16.7 

2.1* 

0.63 

7.0 

- 2 . 1 * 

- 2 . 1 * 













































TABLE IV: Head and tail of angular accelerations. 


«(°) 

pn 

7(“) 

dJAB (rad/s^) 

WflC (rad/s^) 

0 

0 

16.7 

4.2 

-14.1 






90 

90 

16.7 

-1.3 

61.7 


IV. DISCUSSION AND CONCLUSIONS 


The relations needed to simulate the motion of an articulated mechanism were presented and the kinematics was 
analyzed numerically. Results for initial and final positions (a = 0° and 90° respectively) are shown in Tab. |lj 
Tab. |nj Tab. Ill and Tab. |IV| At the initial position (a = 0°) both Vb and ujab have the same numerical values of 
the disk’s rotation rate (cells with asterisks). Paper and pencil solutions give instantaneous values while numerical 
experiments capture the motion’s evolution. Hence, if only a paper and pencil solution is attempted for a = 0°, one 
may be tempted to believe these parameters always have the same numerical values. The same would be true if one 
attempts a paper and pencil solution for the hnal position {a = 90°), since ojab and ujbc have the same absolute 
values of the constant linear speed of a point at the edge of the disk (cells with asterisks). 

The experiments also demonstrated the inversion of rotation for the linking bar under plane motion for two different 
scenarios: one where the rate of rotation of the disk is constant, and one where it is made to vary linearly with the 
angle of rotation. Analysis such as the one presented here complement routine paper and pencil solutions which alone 
may lead to erroneous conclusions about the motions. A design modification of the original mechanism is shown in 
Fig.m The distance between the fixed centers of rotation is shorter than the original mechanism to allow full disk 
revolution while keeping the linking bar size unchanged, and the disk rotation is now clockwise. The full revolution 
is not possible in the original configuration (Fig. [^. 

Despite the complexities of the rigid bodies’ formulations described in this article, one should keep in mind that 
this is still a relatively simple mechanics problem. All principles applied are exact and the numerical implementations 
only required algebraic expressions. Complications arise if one attempts to formulate and model collisions between 
rigid bodies. Impulse-momentum principles are needed compared to the limited kinematics toolbox. Moreover, 
collision laws are not well understood and numerical implementations need to handle differential equations with 
discontinuities that may lead to energy conservation inconsistencies and other difficulties m, 0,0)- 



FIG. 4: Articulated mechanism design modification that allows for full disk revolution. 
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